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ABSTRACT 



Aims. The self-accelerating braneworld model (DGP) seems to provide a simple alternative to the the standard ACDM cosmology 
to explain the current cosmic acceleration, which is strongly indicated by measurements of Type la supemovae, as well as other 
concordant observations. 

Methods. In this work, we investigate observational constraints on this scenario from gravitational lensing statistics using the Cosmic 
Lens All-Sky Survey (CLASS) lensing sample. 

Results. We show that a large parameter space of the DGP model is in good agreement with this radio source gravitational lensing 
sample. 

Conclusions. In the flat case, fl K = 0, the likelihood is maximized, £. = for Q M = 0.30*2 }'. If we relax the prior on Q K > the 

likelihood peaks at {fi M ,fl rt ) ^ (0.29,0.12), just slightly in the region of open models. However the confidence contours are pretty 
elongated so that we can not discard either close or flat or open models. 

Key words, cosmological parameters — cosmology: theory — gravitational lensing — quasars: general 



1. Introduction 

The accelerating expansion of our universe was first discov- 
ered by the measurements of distant Type la supernovae (SNe 
la; Riess et al. 1998; Perlmutter et al. 1999), and was con- 
firmed by the observations of the cosmic microwave background 
anisotropies (WMAP: Bennett et al. 2003) and the large scale 
structure in the distribution of galaxies (SDSS: Tegmark et al. 
2004a,b). By assuming General Relativity, a dark energy com- 
ponent has been invoked as the most feasible mechanism for 
the acceleration. However, although fundamental for our under- 
standing of the Universe, its nature (as well as the nature of the 
dark matter) remains a completely open question nowadays. 

Among several alternatives to dark energy, the models that 
make use of the very ideas of branes and extra dimensions 
to obtain an accelerating universe are particularly interesting 
(Randall and Sundrum 1999a,b). The general principle behind 
such models is that our 4-dimensional universe would be a 
brane embedded into a higher dimensional spacetime bulk on 
which gravity can propagate. One famous brane world model 
is proposed by Dvali, Gabadadze and Porrati (2000), which is 
widely referred to as DGP model. This scenario describes a self- 
accelerating 5-dimensional brane world model with a noncom- 
pact, infinite-volume extra dimension in which the dynamics of 
gravitational interaction is governed by a competition between 
a 4-dimensional Ricci scalar term, induced on the brane, and an 
ordinary 5-dimensional Einstein-Hilbert action. For scales below 
a crossover radius r c (where the induced 4-dimensional Ricci 
scalar dominates), the gravitational force experienced by two 
punctual sources is the usual 4-dimensional 1 /r 2 force whereas 
for distance scales larger than r c the gravitational force follows 



the 5-dimensional 1/r 3 behavior. The Friedmann equation is 
modified as follows 



Ok(1 + if + ( + ^Arc +0 M (1+Z) 3 ) 



(1) 
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where H is the Hubble parameter as a function of redshift z (Ho is 
its value at the present), Q.^, ^r t and £2m represent the fractional 
contribution of curvature, the bulk-induced term and the matter 
(both baryonic and nonbaryonic), respectively. Q It is defined as 
Q, c = l/4r 2 //, 2 . From Eq.(l), the DGP model is a testable sce- 
nario with the same number parameters as the standard ACDM 
model. 

The advantages of the DGP model has triggered a wave of 
interests aiming to constrain its model parameters using vari- 
ous cosmological observations, such as the magnitude-redshift 
relation of supernovae of type la (Avelino and Martins 2002; 
Deffayet et al. 2002; Zhu and Alcaniz 2005; Maartens and 
Majerotto 2006; Barger et al. 2007; Movahed et al. 2007), the 
cosmic microwave background shift parameter from WMAP 
and the baryon acoustic oscillation peak from SDSS (Guo et al. 
2006; Lazkoz et al. 2006; Rydbeck et al. 2007; He et al. 2007), 
the angular size - redshift data of compact radio sources (Alcaniz 
2002), the age measurements of high-z objects (Alcaniz, Jain and 
Dev 2002), the lookback time to galaxy clusters (Pires, Zhu and 
Alcaniz 2006), the optical gravitational lensing surveys (Jain et 
al. 2002), the observed Hubble parameter H(z) data (Wan, Yi and 
Zhang) and the large scale structures (Multamaki et al. 2003; Lue 
et al. 2004; Koyama and Maartens 2006; Song et al. 2007) (For 
a recent review on the DGP phenomenology, see Lue 2006). 

In this paper, we shall consider the observational constraints 
on the parameters of the DGP model arising from the Cosmic 
Lens All-Sky Survey (CLASS) lensing sample. Our results are 



2 



Zong-Hong Zhu and Mauro Sereno: Testing the DGP model with lensing statistics 



in agreement with other recent analyses, providing a comple- 
mentary test to the DGP model. 

Gravitational lensing has been becoming a useful tool for 
modern astrophysics. It provides cosmological tests in several 
ways, such as gravitational lensing statistics (Kochanek 1996; 
Zhu 1998; Cooray & Huterer 1999; Chiba and Yoshii 1999; 
Chae et al. 2002; Sereno 2005), weak lensing surveys (Benabed 
and Bernardeau 2001), Einstein rings in galaxy-quasar sys- 
tems (Yamamoto & Futamase 2001), clusters of galaxies acting 
as lenses on background high redshift galaxies (Sereno 2002; 
Sereno and Longo 2004; Sereno 2007), and gravitational lens 
time delay measurements (Schechter 2004). Results from tech- 
niques based on gravitational lensing are complementary to 
other methods and can provide restrictive limits on the accelera- 
tion mechanism. The aim of the current paper is to check the va- 
lidity of the DGP model with radio-selected gravitational lensing 
statistics. We adopt the Cosmic Lens All-Sky Survey (CLASS) 
statistical data which consists of 8958 radio sources out of which 
13 sources are multiply imaged (Browne et al. 2003; Chae et al. 
2002). We work only with those multiply imaged sources whose 
image-splittings are known to be caused by single early type 
galaxies, which reduces the total number of lenses to 10. We 
show that a large parameter space of the DGP model is in good 
agreement with this radio source gravitational lensing sample. 
The maximum likelihood happens at {Qm,^} - {0.29,0.12}, 
just slightly in the region of open models. 

The paper is organized as follows. In Section 2, the ba- 
sics of gravitational lensing statistics is introduced. Properties 
of the CLASS sample and its statistical analysis are illustrated 
in Section 3. Finally, we present our conclusions and discussion 
in Section 4. 

2. Basics of gravitational lensing statistics 

A realistic statistics of gravitational lenses can be performed 
based on simple assumptions (Kochanek 1996; Chae 2003; 
Sereno 2005; and references therein). The standard approach is 
based on the observed number count of galaxies and on the sim- 
ple singular isothermal sphere (SIS) model for lens galaxies. 

The differential probability of a background source to be 
lensed by a background galaxy with velocity dispersion between 
cr and cr + dcr and in the redshift interval from Zd to Zd + dzd is 



d 2 T driQ cdt 
, , = — (zd,o-)s cr (o-)— , 
dzddcr dcr dzd 



(2) 



where s ct is the cross section for lensing event and is the dif- 
ferential number density of the lens population. For a conserved 
comoving number density of lenses, nc(z) = «o(l + z) 3 - 

The lens distribution can be modeled by a modified 
Schechter function of the form (Sheth et al. 2003) 



driQ I ( 



i 



r(a/j3)cr' 



(3) 



where a is the faint-end slope, B the high-velocity cut-off and n„ 
and cr, are the characteristic number density and velocity dis- 
persion, respectively. Early-type or late-type populations con- 
tribute to the lensing statistics in different ways and type-specific 
galaxy distributions are required. As a conservative approach, 
we do not consider lensing by spiral galaxies. In fact the de- 
scription of the late-type galaxy population is plagued by large 
uncertainties and they contribute no more than 20-30% of the 
total lensing optical depth. A proper modeling of the distribu- 
tion of the lensing galaxies is central in lensing statistics. In 



our analysis we will use the results of Choi et al. (2007) who 
analyzed data from the the SDSS Data Release 5 to derive the 
velocity dispersion distribution function of early-type galaxies. 
They found n, = 8.0xl0~ 3 /i 3 Mpc~ 3 , where h is Hq in units of 
100 km s _1 Mpc-\ cr, = 144+5 km s _1 , a = 2.49 + 0.10, and 
B = 2.29 + 0.07. 

Early-type galaxies can be well approximated as singu- 
lar isothermal spheres. As shown in Maoz & Rix (1993) and 
Kochanek (1996), radial mass distribution, ellipticity and core 
radius of the lens galaxy are unimportant in altering the cos- 
mological limits. Assuming a flat model of universe, a typical 
axial ratio of 0.5 in a mixed population of oblate and prolate 
spheroids would induce a shift of ~ 0.04 in the estimation of Hm 
(Mitchell et al. 2005), well below statistical uncertainties. Since 
departures from spherical symmetry induce a relatively small ef- 
fect on lens statistics and the distribution of mass ellipticities is 
highly uncertain, spherically symmetric models supply a viable 
approximation. The cross section of a SIS is 



16tt 



(4) 



where Dd, £>ds and D s are the angular diameter distances be- 
tween the observer and the deflector, the deflector and the source 
and the observer and the source, respectively. The two multiple 
images will form at an angular separation 



A<9 = 8 



\c) D s 



(5) 



which relates the image separation to the velocity dispersion of 
the lens galaxy. The total optical depth for multiple imaging 
of a compact source, t, the probability that a SIS forms mul- 
tiple images of a background source with angular separation 
AO, dr/dAO, and the probability of lensing by a deflector at Zd, 
dr/dzd, can be obtained by integrating the differential probabil- 
ity in Eq. (fJJ. 

Lensing probabilities must be corrected for the magnifica- 
tion bias B, i.e. the tendency of gravitationally lensed sources 
to be preferentially included in flux-limited samples due to their 
increased apparent brightness (Turner 1990; Fukugita & Turner 
1991 ; Fukugita et al. 1992; Kochanek 1993). The bias factor for 
a source at redshift z s with apparent magnitude m is given by 



B(m,z,M ) = 

\ dm 



(6) 



r + °° dN s 

Jm dm 



(m + 2.5 log M, z)P(M)dM, 



A/o being the minimum magnification of a multiply imaged 
source, with value M = 2; P(M)dM = 2M\M^dM is the 
probability that a multiple image-lensing event causes a total 
flux increase by a factor M (Kochanek 1993). The function 
dN s /dm is the differential source number count in magnitude 
bins dm. Furthermore, since observations have finite resolution 
and dynamic range, lens discovery rates are affected by the abil- 
ity to resolve multiple source images (Kochanek 1993). Lensing 
probabilities must then account for the resolution limit of the sur- 
vey. For the SIS model, selection effects can be characterized 
by the maximum magnitude difference that can be detected 
for two images separated by A6, Am(A6), which determines a 
minimum total magnification M t = Mo(f + l)/(/ - 1), where 
2.5 log / = Am (Kochanek 1993). 
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Finally, the likelihood function can be written as (Kochanek 
1993; Chae et al. 2002) 



■c=Y\a-pdY\pi,j, 

7=1 



(7) 



;=i 



where A^l is the number of multiple-imaged sources and Ny is 
the number of unlensed sources, pi is the suitable probability 
accounting for the whole of the data available for each lens sys- 
tem, i.e. the lens redshift and/or the image separation (Chae et al. 
2002; Mitchell et al. 2005). Probabilities are corrected for bias 
and selection effects. 

Since r <k 1 the likelihood can be approximated as (Mitchell 
et al. 2004) 



X ^ exp 



N z (zs)p(Zs)dz s 



(8) 



where N z (z s ) is the redshift distribution of the sources. We use a 
uniform distribution for the priors on the cosmological parame- 
ters, so that, apart from an overall normalization factor, the like- 
lihood can be identified with the posterior probability. 

3. Data analysis 

In this section, we discuss the radio-survey used for our lensing 
statistics and present the constraints on the parameters of the 
DGP model. 

3.1. Data set 

The most reliable data set suitable for statistical analysis is 
provided by a sample of 8958 flat-spectrum radio sources 
with 13 lenses by the Cosmic Lens All-Sky Survey (CLASS; 
Browne et al. 2003; Myers et al. 2003). Data of interest are 
listed in table 1 of (Chae 2005). We limit our analysis to 
the early-type lens galaxies. Ten systems in the CLASS sam- 
ple (0445+123, 0631+519, 0712+472, 1152+199, 1359+154, 
1422+231, 1608+656, 1933+503, 2114+022 and 2319+051) 
can be assumed to be early-type lenses (Chae 2005). We do not 
consider the information on the image separation in 1359+154, 
1608+656 and 2114+022 whose splittings are strongly affected 
by galaxy companions very close to the main lens. 

The final CLASS statistical sample has been selected such 
that, for doubly imaged systems, the flux ratio is < 10 and it 
is independent of the angular separation. According to the se- 
lection criteria, the compact radio-core images have separations 
greater than = 0.3 arcseconds. The probabilities that en- 

ter the likelihood must be then considered as the probabilities 
of producing image systems with separations > A# m ; n . Taking 
into account the CLASS observational selection function, Chae 
(2007) found a magnification bias of B 3.36 for the SIS. 

Redshift measurements are only available for a restricted 
CLASS subsample. Following Sereno (2005), we model the red- 
shift distribution N z (z s ) of the sources with a kernel empirical 
estimator. For the unmeasured lensed source redshifts, we set z s 
to the mean redshift of the sources lensed by early-type galaxies 
with measured redshift, (z s )i en sed = 2.2. 

3.2. Statistical analysis 

Let us now perform a statistical analysis of the data sample. As a 
first step, we fix the nuisance galactic parameters to their central 
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Fig. 1. Normalized likelihood, X/X n 
a flat geometry, £2k = 0. 
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Fig. 2. Normalized likelihood, X/Xmax, in the QM-£2r c plane. 
The dot shows the best fit model and the contours denote the 
68.3%, 95.4% and 99.7 confidence limits for two parameters. 
The dashed line represents the locus of flat models of universe 
(£2k = 0); bouncing models in the upper-left shaded region do 
not have big bang. 



values. We will consider the related uncertainty later. In the flat 
case, £2k = 0, the likelihood is maximized, X = Xmax, for Q.m = 
0.30*9 Jp see Fig.[T] Uncertainties denote the statistical 68.3% 
confidence limit for one parameter, determined by X/Xmax = 
exp(-l/2). 

Even if we relax the prior on Qk, the likelihood peaks 
for nearly flat models. In fact, the likelihood is maximum for 
{Om^J — {0.29,0.12}, just slightly in the region of open 
models, see Fig. [2] The three contours in the figure correspond 
to the 68.3%, 95.4% and 99.7% confidence limits for two pa- 
rameters, namely X/Xmax = exp(-2.30/2) , exp(-6.17/2) and 
exp(- 11.8/2), respectively. However, contours are pretty elon- 
gated so that we can not discard either close or flat or open mod- 
els. 

Uncertainties in the redshift distribution of the sources can 
induce additional errors in the estimates of the cosmological pa- 
rameters. A source of error is the finite sample size of the sam- 
ple of measured source redshifts (only 27 source redshifts are 
known), which induces an error in the estimated redshift distri- 
bution. From a bootstrap resampling procedure, it can be cre- 
ated a set of simulated distributions which is then used to cre- 
ate a new kernel estimator for the redshift distribution. It can be 
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shown that the finite size induces a dispersion of ~ 0.08 on Qm 
(Sereno 2005). On the other hand, the cosmological constraints 
are nearly insensitive of the functional form used when model- 
ing the redshift distribution. Conclusions are really unaffected if 
a Gaussian distribution is used instead of the kernel estimator. 
Finally, results change in a very negligible way if we use differ- 
ent values of z s for the lensed sources with unknown redshift. 

The main uncertainty in the estimation of cosmological pa- 
rameters comes from errors in the assumed parameters of the ve- 
locity dispersion distribution function which describes the lens 
population. In order to estimate such source of error, we simu- 
lated a sample of 100 sets of galactic parameters by extraction 
from normal distributions centered on the best estimates of each 
parameter and with standard deviation given by the associated 
uncertainty. The likelihood analysis was then repeated for each 
set of galactic parameters. Assuming flat cosmological models, 
the resulting distribution of the maximum likelihood estimates 
has a scatter of ~ 0.09, which gives a similar uncertainty in the 
determination of Qm- 

Finally a theoretically important systematic uncertainty is 
due to the the effect of small-scale inhomogeneities on large- 
scale observations. Matter distribution is locally inhomogeneous 
and affects light propagation and the related cosmological dis- 
tances (Sereno et al. 2001; Sereno, Piedipalumbo and Sazhin 
2002; and references therein). However, being the universe glob- 
ally homogeneous, the effect on the total lensing statistics is 
small (Covone et al. 2005). 

4. Conclusion and discussion 

Since the discovery of the accelerating expansion of the uni- 
verse, in addition to the standard ACDM cosmological model, 
a huge number of scenarios have been proposed to be the 
acceleration mechanism (for a recent review, see: Sahni and 
Starobinsky 2000; Padmanabhan 2003; Lima 2004; Copeland, 
Sami and Tsujikawa 2006; Alcaniz 2006). Examples include the 
so-called "X-matter" (Turner and White 1997; Zhu, Fujimoto 
and Tatsumi 2001; Alcaniz, Lima and Cunha 2003; Dai, Liang 
and Xu 2004; Rupetti et al. 2007; Wang, Dai and Zhu 2007), 
a decaying vacuum energy density or a time varying A-term 
(Ozer and Taha 1987; Vishwakarma 2001), an evolving scalar 
field, dubbed quintessence (Ratra and Peebles 1988; Caldwell 
et al. 1998; Wang and Lovelace 2001; Gong 2002; Chen and 
Ratra 2004; Choudhury and Padmanabhan 2005; Ichikawa et 
al. 2006), the phantom energy, in which the sum of the pres- 
sure and energy density is negative (Caldwell 2002; Dabrowski 
et al. 2003; Wang, Gong and Su 2004; Wu and Yu 2005, 2006; 
Chang et al. 2007), the Chaplygin gas (Kamenshchik et al. 2001 ; 
Bento et al. 2002; Alam et al. 2003; Alcaniz, Jain and Dev 
2003; Dev, Alcaniz and Jain 2003; Silva and Bertolami 2003; 
Makler et al. 2003; Zhu 2004; Zhang and Zhu 2006), the quin- 
tom model (Feng, Wang and Zhang 2005; Guo et al. 2005; Zhao 
et al. 2005; Xia et al. 2006; Wei and Zhang 2007), the holo- 
graphic dark energy (Li 2004; Zhang and Wu 2005; Chang, Wu 
and Zhang 2006), the Cardassion model (Freese and Lewis 2002; 
Zhu and Fujimoto 2002, 2003; Sen and Sen 2003; Wang et al. 
2003; Gong and Duan 2004a,b; Wang 2005; Bento et al. 2006; 
Reboul and Cordoni 2006; Yi and Zhang 2007) and the Casimir 
force (Szydlowski and Godlowski 2007; Godlowski et al. 2007). 
All these acceleration mechanisms should be tested with various 
astronomical observations. 

In this paper, we have focused our attention on the DGP 
model. We have analyzed this scenario by using the Cosmic 
Lens All-Sky Survey sample (Browne et al. 2003; Myers et al. 



2003) to obtain the 68.3%, 95.4% and 99.7% confidence re- 
gions on its parameters. It is shown that a large parameter space 
of the DGP model is consistent with this radio source gravita- 
tional lensing sample. In the flat case, Qk = 0, the likelihood 
is maximized, £. = -£max, for £2m = 0.30^' [j*. If we relax the 
prior on Q.^, the likelihood peaks at |£2m,£\ ) - {0.29,0.12}, 
just slightly in the region of open models. The obtained con- 
fidence regions of Figure 2 are also in good agreement with 
the results from analyzing data of type la supernovae (Zhu and 
Alcaniz 2005), which implies that gravitational lensing statis- 
tics provides an independent and complementary constraint on 
the DGP model. However, similar to the case of type la su- 
pernovae, the confidence contours are pretty elongated so that 
we can not discard either close or flat or open models by only 
using the CLASS sample. Using the gold sample of type la 
supernovae (SNela), the first year data from the Supernova 
Legacy Survey (SNLS) and the baryon acoustic oscillation 
(BAO) peak found in the Sloan Digital Sky Survey (SDSS), 
Guo et al (2006) obtained, at 99.73% confidence level, Q. m = 
0.270!™!* and Q rc = 0.216+°;°}* (hence a spatially closed uni- 
verse with Qi = -0.350^}yg3), which seems to be in contra- 
diction with the most recent WMAP results indicating a flat 
universe. Based on this result, the authors also estimated the 
transition redshift (at which the universe switches from de- 
celeration to acceleration) to be 0.70 < z q =o < 1-01, at 2<x 
confidence level. Therefore, the method of combining obser- 
vational data provides much more stringent constraint on 
the DGP model than any single data. It is naturally hopeful 
that, with either future larger gravitational lensing samples or 
a joint investigation with other astronomical observations, one 
could obtain a more stringent constraint on the DGP model pa- 
rameters. 
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